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We report the first experimental observation of noise-free 
stochastic resonance by utilizing the intrinsic chaotic dynam- 
ics of the system. To this end we have investigated the effect 
of an external periodic modulation on intermittent signals ob- 
served by high power ferromagnetic resonance in yttrium iron 
garnet spheres. Both the signal-to-noise ratio and the resi- 
dence time distributions show the characteristic features of 
stochastic resonance. The phenomena can be explained by 
means of a one-dimensional intermittent map. We present 
analytical results as well as computer simulations. 

PACS numbers: 05.45.+b, 05.40.+j, 75.30.Ds 



I. INTRODUCTION 

Stochastic resonance, invented 15 years ago as a model 
for geophysical dynamics 0], has meanwhile found its 
way into such diverse fields like physics, meteorology, 
chemistry, and biology [^| |J. The rising interest in this 
field stems from the counterintuitive effect that a peri- 
odic signal component can be amplified by a stochastic 
force. 

The basic mechanism is usually explained in terms 
of the famous Kramers problem S, i.e. the overdamped 
motion of a particle in a symmetric double-well potential 
subjected to noise, which is supplemented by a time pe- 
riodic forcing. The noise causes "incoherent tunnelling" 
between the two wells with an exponentially decreasing 
distribution of the respective residence times. The pe- 
riodic forcing leads to enhanced transitions on certain 
time scales and, therefore, to a periodic signal compo- 
nent. It is the prominent feature of stochastic resonance 
that the signal-to-noise ratio does not decrease with in- 
creasing noise amplitude, but attains a maximum at a 
certain noise strength. A second characteristic property 
shows up in the distribution of residence times, where 
the periodic forcing leads to maxima at odd multiples 
of half the driving period (cf. ||). Of course, these sig- 
natures of stochastic resonance are not confined to this 
special model, but occur in more general bi- and monos- 
table systems and for different types of noise HQ. 



It was suggested by Anishchenko et al. Q| based on 
the numerical analysis of a simple map, that quite sim- 
ilar phenomena can be caused by the internal chaotic 
dynamics instead of an external noise. In that case the 
intermittent hopping between different chaotic repellers 
in conjunction with an additional periodic forcing can be 
used to amplify the periodic signal component in close 
analogy to conventional stochastic resonance. This math- 
ematical model system motivated us to look for noise- 
free stochastic resonance in real experimental systems. 

In this letter we report on its realization in chaotic 
spin-wave dynamics. By taking advantage of the known 
bifurcation scenarios beyond the instability threshold 
||f§, we have investigated the chaotic intermittent dy- 
namics subjected to an additional periodic forcing. 

II. EXPERIMENTS AND RESULTS 

High power ferromagnetic resonance experiments were 
performed on spheres of yttrium iron garnet (YIG), 
which is well established as a "prototype nonlinear 
ferromagnet" and has extensively been studied for 
decades [[10-12 . The sample was placed in a bimodal 
transmission-type microwave cavity and excited by a mi- 
crowave field of 9.25 GHz, applied perpendicularly to the 
static magnetic field H. The parametric excitation of 
spin waves was observed in subsidiary absorption [ jl3| . 
The transmitted and rectified signal was recorded with 
a digital transient recorder and analysed on a PC. Ac- 
cessible system parameters are the static magnetic field 
and the microwave power, which is proportional to the 
squared amplitude of the pumping field. Increasing the 
pumping power above the first-order Suhl threshold fTs] ] 
we observed auto-oscillations and a variety of routes into 
chaos: period doubling, quasi-periodicity, and different 
types of intermittency 1 1 1 15] . Our measurements were 
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performed in a parameter region of type-Ill intermit- 
tency, where the system jumps between period 2 be- 
haviour (regular phase) and chaotic behaviour (chaotic 
phase). The intermittency scenario can be run through 
in two ways, either by varying the microwave input power 
P or the static magnetic field H |J. 

Here, the regular and chaotic phases take the role 
of the two states in conventional models of stochastic 
resonance, and the intrinsic chaotic dynamics acts as a 
"stochastic" forcing. Accordingly, the mean lengths of 
the two phases can be changed on variation of the two 
external control parameters P and H. In this sense P 
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and H correspond to the noise strength in conventional 
stochastic resonance experiments. To obtain the reso- 
nance effect one has to apply an additional periodic 
forcing, which in our case consists of an amplitude mod- 
ulation of the microwave power. The two additional pa- 
rameters, modulation frequency and amplitude, have to 
be adjusted properly, i.e. the modulation amplitude has 
to be chosen small enough to ensure that the system stays 
inside the intermittency regime, but large enough to ob- 
serve a periodic component in the transmitted signal. 
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FIG. 1. Mean lengths of the chaotic (A) and the regular 
(V) phases vs. modulation frequency. 
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FIG. 2. (a) Probability distribution of the chaotic lengths 
for v — Vmax- The full line indicates the exponential back- 
ground, (b) Difference between the total probability and the 
exponential background taken from (a). 

The influence of modulation frequency shows up in the 
duration of regular and chaotic phases (cf. fig.|l|). Wc 
found that the mean length of the regular state is not 



affected by the modulation. The mean chaotic length 
shows a distinct maximum at v max — 1kHz, which means 
that the reinjection from the chaotic to the regular state 
becomes less probable. At v max the mean return time, 
i.e. the sum of the two lengths, is of the order of the mod- 
ulation period, which is a prerequisite for the occurrence 
of stochastic resonance. 

Signatures of stochastic resonance could be observed 
in the distribution of the durations of the chaotic phase, 
i.e. the residence time distribution. Without modula- 
tion one would obtain an exponential decay, caused by 
the uniform reinjection from the chaotic to the regular 
state. With modulation there appears a structure on top 
of this background, which represents a typical feature of 
stochastic resonance [|| . The distribution at v max is plot- 
ted in fig. ||(a). The exponential background resulting 
from the intermittent dynamics is larger than in conven- 
tional stochastic resonance. Figure ||(b) shows the dis- 
tribution after having subtracted this background. The 
distance between the peaks is equal to the modulation 
period T m , as expected from theory. Since our system is 
strongly asymmetric the first peak is not located at ex- 
actly half the period but shifted towards a higher value. 
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FIG. 3. Signal-to-noise ratio vs. the bifurcation parame- 
ters ei = P c — P (a) and Ei = H c — H (b) for various modu- 
lation frequencies ( □ 1.0kHz, o 5.0kHz, A 10.0kHz ). 

Stochastic resonance was originally defined by an in- 
crease of the signal-to-noise ratio with increasing noise 
strength. In our experiments we have two possibilities to 
change the internal " deterministic noise level" . At the 
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bifurcation point (H = H c , P = P c ) the regular state is 
marginally stable. With increasing bifurcation parame- 
ters E\ = P c — P and E2 = H c — H Q the switching rate 
between regular and chaotic states increases and hence 
the internal noise level, too. If we take, as usual, the am- 
plitude ratio between peak and background of the Fourier 
spectrum as a definition of the signal-to-noise ratio, we 
observe a maximum on variation of the static magnetic 
field or the microwave power. Figure ^ shows the results 
for several modulation frequencies on a logarithmic scale. 
Distinct resonance phenomena can be seen in both the 
field and power dependences. 



ent values 73 and 73 every (n p /2) th time step. Let pi (n) 
denote the probability density for a trajectory hitting the 
interval /; at time n. It is periodic in time with period 
n p but attains a constant value on each interval, since 
the model maps intervals on intervals (Markov map). 
Furthermore let Sk(n) C I denote those phase space 
points which stay exactly k time steps in the chaotic re- 
gion I\ U I2 U I3 if the initial phase of the modulation 
is n, and denote the size of this set by X(Sk(n)). Then 
pa(n)\(Sk{n)) gives the probability that at phase n of the 
modulation a chaotic burst of length k occurs. Hence by 
averaging over the initial phase we obtain the desired 
distribution of residence times 



III. THEORETICAL MODEL AND 
SIMULATIONS 

Far above the Suhl threshold an approach from first 
principles has turned out to be extremely complicated 
|TT|| . Thus we adopt a phenomenological description in 
terms of simple mathematical model systems, which have 
proven to be fruitful in diverse problems of nonlinear dy- 
namics. In order to develop a theoretical explanation for 
the residence time distribution we propose a time depen- 
dent one-dimensional map which can be treated analyt- 
ically. 



1 



T(x) 









'3\\ 




/ I I I 



I. 







I, 



FIG. 4. Piecewise linear map x n +i = T(x n )- The broken 
line in I3 and the double arrow show the time dependent 
modulation. The shaded box in the upper right corner indi- 
cates the region of the chaotic repeller. 

This map generates a regular phase via an inverse 
pitchfork bifurcation and a chaotic phase via a chaotic 
repeller (cf. fig|I|). We adopt a piecewise linear model 
with slopes 7/ on the intervals // to keep the analysis en- 
tirely analytical ]T^ | . A periodic modulation with period 
n p is included via the time dependence of the slope in the 
right hand interval I3. For simplicity we confine ourself 
to the case where the slope 73 changes between two diffcr- 
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n pPo X(I ) j- 



po{n)\(S k (n)) 



(1) 
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where po denotes the average of po(n) over one period. 
The structure of the distribution is determined by the 
time dependence of the density as well as by the time 
dependence of Sk (n) . The size of S k (n) can be expressed 
by the (static) escape rates of the chaotic repeller |j[5| . 
Using the abbreviation 



PA := l/M + l/M + VI-tfl 
_ . if (k + n)|mod n p < n p /2 

p B := 1/I71I + 1/I72I + 1/I73 3 ! 
if (k + n)|mod n p n p /2 

a geometrical consideration yields 
fc-i 

X(S k (n)) = A(7 ) l[pi(n) ■ (1 -j»*(n)) 



(2) 



(3) 



The density po{n) can in principle be evaluated, too. But 
to leading order in the modulation amplitude it remains 
time independent po(n) — po- This implies that the prop- 
erties of the regular phase are not affected by the mod- 
ulation, which is in accordance with our experimental 
observation. With eq.(^) the distribution |l]) reads 



-Ell^w-d 



Pk(n)) 



(4) 



Its structure is entirely determined by the time depen- 
dence of the escape rate and does not depend on details 
of the intermittency mechanism. Expression (^) is evalu- 
ated by combinatorics. In the case of large period n p 3> 1 
the distribution becomes quasi-continuous 



N k ~ exp(-crfc/n p ) | 



f-(x), 0^x<l/2 
f+(x), 1/2 ^x<l 



(5) 



where x = fc/n p |mod 1 denotes the fractional part of 
the chaotic length with respect to the modulation pe- 
riod. Expression (||) reflects an exponential decay of the 
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distribution with rate a := n p [l — (pa + Pb)/2] and a 
superstructure of period n p (cf. figj^) , which is explicitly 
given by 

/±(l/2 ± z) := 2zo-cosh[<5(l/2 - z)\ 
±2^sinh[(5(l/2-^)] + 2cr/(5sinh[(5(l/2-z)] . (6) 

Here 5 := (j>b — PA)n p /2 represents the amplitude of the 
modulation. 

In addition to our analytical results we have performed 
computer simulations using a one-dimensional map like 
fig. [| but with the left-hand part being replaced by a 
cubic function. The time dependence was introduced 
by modulation of the right-hand maximum. From the 
time series we evaluated the distribution of the residence 
times which are shown in fig. [5]. Like for the analyti- 
cal results we found an exponential decay modulated pe- 
riodically with the period of the forcing. In this case, 
like for the ferromagnetic resonance experiments, a large 
exponential background occurs which is caused by the 
internal chaotic dynamics, i.e. the "deterministic noise 
level" . The coincidence is even more convincing when 
keeping in mind that the analytical prediction and the 
simulation have been obtained from different maps and 
different modulation mechanisms. 



the signal-to-noise ratio and the residence time distri- 
butions show exactly the same characteristic behaviour 
as in conventional stochastic resonance. We were able 
to explain these features in terms of a one-dimensional 
intermittent map. An analytical expression for the res- 
idence time distribution was derived and compared to 
results from computer simulations. Both the quan- 
titative agreement between theoretical calculations and 
simulations, and the qualitative coincidence with the ex- 
perimental result demonstrates the universality of these 
features. We expect that noise-free stochastic resonance 
will be realized in many other physical or technical ap- 
plications in the near future. 
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FIG. 5. Distribution of chaotic lengths. Analytical result 
(solid line) and simulation. 

IV. CONCLUSIONS 

We have reported on the experimental realization of 
stochastic resonance in ferromagnetic resonance experi- 
ments beyond the Suhl threshold. In contrast to conven- 
tional stochastic resonance we used the interplay between 
the intrinsic chaotic dynamics in an intermittent parame- 
ter region and an external periodic modulation for the re- 
alization of noise-free stochastic resonance. Nevertheless 
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